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Long-ranged electrostatic interactions in electrolytes modify their contact angles on charged sub- 
strates in a scale and geometry dependent manner. For angles measured at scales smaller than 
the typical Debye screening length, the wetting geometry near the contact line must be explicitly 
considered. Using variational and asymptotic methods, we derive new transcendental equations for 
the contact angle that depend on the electrostatic potential only at the three phase contact line. 
Analytic expressions are found in certain limits and compared with predictions for contact angles 
measured with lower resolution. An estimate for electrostatic contributions to line tension is also 
given. 

PACS numbers: 47.10.-|-g, 68.08.-p, 68.08.Bc 



Modern microfluidic and patterning applications call 
for directed fluid flow and wetting on treated surfaces 
often exhibiting complex surface chemistry ^, ||. 
Mechanisms of differential wetting of small droplets of 
electrolytes have also been exploited as electrically ac- 
tivated switches and micropumps In such systems, 
the effects of surface ionization have been found to be 
important |, |, @, |. 

Consequently, one often considers two immiscible flu- 
ids of dielectric constants eo,ei that partially wets an 
ionizable, rigid substrate, as shown in Figs. 1. The 
two liquids, depending on their individual pKa's, will 
differentially hydrolyze/ionize substrates such as glass. 
Ionizable surfactants may also be adsorbed, imparting 
a fixed, relatively uniform surface charge a at the in- 
terfaces. Microscopically, surface tensions 7 arise from 
mismatches in short-ranged molecular interactions {e.g. 
van der Waal's) among the various species. Electrostatic 
double layers also contribute to surface energies. The 
application of the classical Young-Duprc (Y-D) equation 
1^, 7giC0sa* ~ — Yis, with electrostatically modi- 
fied surface energies {e.g. 7^^ = 71s -|~ ^(Tiipi, where cti 
and ipi, are surface charges and potentials in, say, liq- 
uid 1 far from the contact point P) is accurate provided 
the apparent contact angle a* is measured at a point 
P* outside the range of the electric double layers. 
However, the screening length k~^, although typically 
smaller than the resolution of optical goniometry mea- 
surements, can be within the resolution (nanometers) of 
emerging angle measurement techniques using AFM M. 
Even at lower resolutions, finite-sized double-layer effects 
may be relevant for contact angle measurements. For ex- 
ample, hydrolysis of pure water gives k^^ ^ 1/im, while 
in organic mixtures, with fewer mobile ions, the screen- 
ing length can be even longer If the contact angle is 
measured at P, within the ionic double layers, the simple 
Y-D equation is not appropriate. 

Previous theories that consider surface energy mod- 
ifications 1^, 1^, 1^, ^, |ri|, |l2|, |l3j, have either assumed 
microscopic-ranged interactions or infinite-system sur- 
face free energy changes {i.e. 7 — > 7 -t- ^crcp). Effects 
of the wedge-like geometry on the intrinsic electrostatics 
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FIG. 1: (a) Two liquids with dielectric constants eo, ei wets a 
charged substrate, (b) The substrate is ionized and acquires 
fixed charge densities of (Tq, cri. Charged surfactants can also 
impart a surface charge am at the fluid-fluid interface, (c) In 
order to satisfy boundary conditions, the different screening 
layers deform near the contact point P. 



near the contact line have not been considered. In this 
letter, we derive formulae for the angle^ at the true three 
phase contact line in the presence of long-ranged electro- 
static interactions. Rich features arise in this most simple 
and classic problem when geometry is self-consistently 



^ Measured within both double layers, but outside the length scale 
of other more microcscopic interactions. 
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incorporated. We propose using weak ionic solutions 
as a system with precisely controllable electrostatics for 
studying charge contributions to contact angle and line 
tension forces. Our results are summarized by Eqns. (|^), 
(0), and 

The "mechanical" free energy of an axisymmetric liq- 
uid droplet of footprint radius R (cf. Figs. la,b) in 
contact with a flat, solid substrate is Gmech ~ [jis ~ 
70s) J dr + 7oi f dSoi + Gbody, where 7 are the surface 
tensions sans electrostatic interactions, dr is the in- 
plane surface element, S is the liquid-liquid surface ele- 
ment. The Gbody term may include gravitational energy, 
(pi — /9o)f//2 J h'^(r)dr, and/or Lagrange multipliers to 
e.g. fix droplet volumes of incompressible liquids. The 
electrostatic free energy for a specified surface charge en- 
semble, written as a functional of the local electrostatic 
potential ip is 



Gee = ^ a{S)ip{S)dS - J drdz 



e(r, z) 



|V^|2 + [/[^] 



(1) 

where C/[(^(r)] = eMiC°°(exp(— :^i(p(r)) — 1) is a term 
summing the interactions among mobile charged species 
i, each with valency Vi and bulk concentration c°°. We 
have expressed all energy and length quantities in terms 
of ksT and the Bjerrum length £b = e^/ksT. Variation 
of G'e£[iy9(r), V(^(r); /i(r), V/i(r)] with respect to (fi{r) for 
a fixed droplet height function h{r) yields the Poisson- 
Boltzmann equation with appropriate boundary condi- 
tions. Similarly, variation of Gt = Gmech + Gei with 
respect to the droplet height h{r) determines the com- 
plete shape of the electrolyte droplet via 



7oi9t/i±(r) 



1 



(2) 



where h± is the deformation normal to the tangent t rel- 
ative to a constant slope. Further minimizing the bound- 
ary terms in Gt (which are independent of C/[</5]) with 
respect to the position of the contact point, SGt^z = 
0)/SR^O, yields^ 



701 cos a = (70s - 7is) + (o-o - cri)ip{P; a). (3) 

The first two terms arise from setting the variations 
in the boundary terms of Gmech to zero and reproduces 
the Young-Dupre equation. The last term in (^ is a 
new generalization of the Y-D equation and arises from 
minimizing the boundary terms of Gei. The additional 
term depends only on the jump in the solid surface 
charge and the electrostatic potential (p{P) at the three 



^ The liquid-liquid surface tension 701 may also change due to 
electrocapillarity, but this can be modification measured inde- 
pendently using e.g. pendant drop methods. 



phase contact line P. This potential is found by solv- 
ing the Poisson-Boltzmann equation W'^ip{r) = J7'[(^(r)] 
in the appropriate geometry, subject to boundary condi- 
tions. Therefore, ip{P) will depend parametrically on the 
droplet shape (and hence the contact angle a) near P. 
Equation (^ is exact provided the electrostatic energy is 
given by Gee , and gives an implicit formula for predicting 
the contact angle a. 

In the following, we compute (p{P; a) in the linearized 
limit {U[(p] « «;W2), valid for eif{r)/kBT < 1, by 
solving Kelvin's {aka Debye-Hiickel) equation in each of 
the two fluid domains depicted in Figs. 1: 



Aifj = K^tpjlr, 9) in fluid = 0, 1. 



(4) 
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are as- 



The screening lengths Kj ^ = (^^^ J2i 

sumed much smaller than the dimensions of the droplet. 
Furthermore, neglecting gravity, (^) shows that /i_L(r) 
varies (relative to a perfect wedge) over a length scale 
Lh ^ ^7oi/(k(p^(P)). Provided Lh 3> k^^, the wedge 
is distorted only in the region where 1^9 « 0, sufficiently 
beyond the contact point P to be appreciably influenced 
by electrostatic interactions. Thus, '^{P) will be com- 
puted using a perfect wedge geometry. The boundary 
conditions associated with (W) in 2D wedge domains are 



ei9e<pi(r, 0) = A-Krui 

(po{r,a) - (pi{r,a) = 
eodeifoir, a) - eide(pi{r, a) = Anra^ 



(5) 



The linear problem defined above is related to the classic 
problem of wave scattering from a wedge, which remains 
a substantial mathematical and computational challenge 
p5[ . The problem is best attacked using the Lebedev- 
Kantorovich (LK) integral transform 

/■°° dr 



and its inverse 
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V]{rid) = ^ / f{x,9)Ki^{Kjr)xsh{-nx)dx = C~lfj. 
I" Jo 

(7) 

If the contact angle can be measured sufficiently 
close the contact point (closer than k^^, such as 
in ijl^), the potential ip{Kr P) can be eas- 

ily expressed in terms of its transform f(x,9). We 
will henceforth scale distance as ^ = K^r and rein- 
sert Ko back into the final, quoted results. Asymp- 
totic analysis of the integral representation Kix{S,) = 
(l/2)X!^cxp[-^chi + ixi]dt yields lim^^o Ki^{^) ^ 

{smxt*)/x, where t* = lim^^o cosh~^(l/^) ~ 00. The 
limit limt*^oo(sina:i*)/a; = 'kS{x), implies 



\im^(^,e) = lim \x^f(x,e)] 



(8) 



Wc find analytic expressions for ip{r, 9) in two limiting 
cases that illustrate the full range of behaviors: kq ~ ^ii 
and Ki 3> kq- In the former case of nearly identical 
screening lengths unequal dielectric constants eo 7^ ei 
provide an nontrivial electric field jump across 9 = a. In 
the latter case, we will choose a small screening length 
{e.g., high salt) in region 1 without loss of generality. In 
both cases, the potential can be expanded in the power 

series (pj{r,9) = J2^=o t^^"\^)V''f\'''y^)^ where e is a 
small parameter that depends on the ratio of screening 
lengths and the relevant regime: 

Ki « kq limit - In this limit, = e" and 

e = (ki/kq)^ ~ 1- At each order e" the governing 
equations are A(/3q"'' ~ i^Sq"' for 9 > a, and A(/?^"'' = 
95^"^ + for 9 < a. Upon applying Ci to the 7i*'' 

order Debye-Hiickel equations, 9|/g"'* = fj^^\x, 9) and 
= (2;, 0) = (x, 0) where ^("-i) (x, 0) = 

£i^("-i)r2(l-5„,o). 

Similarly, the transformed boundary condi- 
tions at each order become dgfQ^\x,9 = tt) = 
(2^Vo/eo)sech(7ra;/2)<5„,o, dgf[''\x,9 = 0) = 
— (27r^cri/ei)sech(7rx/2)(5„_o, and eidefi"\x,9 = 
a) — eodefQ^\x,9 = a) = 27r^(JmSech(7ra;/2)i5„,o- Since 
f^\x,9) obeys a homogeneous equation, it is deter- 
mined explicitly. The equations for f^^\x,9) can be 
solved with the appropriate Green function G{x,9,9') = 
— chx6'<cha;(0> — o.)/{xsh.xa),de'G{x,9,9' = 0, a) = 0, 
where 0<^{9'-^) is the smaller (larger) of 9, 9' . We obtained, 
after some algebra, an integral equation for fi(x,9). 
The first iteration of this integral equation yields 

/•oc pa 

ef[^\x,9)^e dx' d9' H{9,9';x,x')f["\x' ,9') 
Jo Jo 

(9) 

where 

H{9,9';x,x') EE 

(ei/4)G(6l, 9'; x)(x^ ~ x'2)a;'sh7ra;'sha;a 
(eishxa -|- eocha;6' tha:(7r — a)) sh^(x + a;')sh^(x — x') 

Upon integrating over 9', and using fi{x,9) w 
f[°\x, 9)+ef[^\x, 9) in (|), we find if{P- a) up to 0(e): 



27r2/Ko 



eia + eo(7r — a) 



(cto + CTl + (Tm) 



-|-£(CTO'5o(a; A) + aiSi{a; A) + (T„iS',„(a; A)) -I- 0(£^) 

(11) 




FIG. 2: The functions Si{a; X = eo/ei) representing the order 
£ — ki/kq — 1 contributions to ip{a; P) (cf. Equation (|ll[)). 



The first (zeroth order) term arises from lim^^^g 
while the 0(e) terms can be expressed as single a- 
dependent integrals Si{a;X). The functions Si are plot- 
ted vs. a in Fig. ^ for various dielectric mismatches 
eo/ei = A. The term that remains when kq = ki and 
that multiplies all higher order terms in e is an effec- 
tive angular average over the dielectrics eojEi- This 
dominant term is qualitatively different from that aris- 
ing from simple surface tension renormalization 7 — > 
7 -|- ^tT(/? = 7-1- 27rcr^/(eK). The higher order terms Si 
are determined with ni relative to kq and thus have 
their major effect when a « tt as more of the volume 
is occupied by electrolyte of inverse screening length ki. 
For larger (smaller) A = eo/cij So{Si) varies more sig- 
nificantly over a smaller range of angles a, reflecting 
the importance of the charged surface on the lower di- 
electric slice. The effects of the liquid-liquid surface 
charge are symmetric with the interchange a it — a, 
eo, kq^ctq ci, Ki, (Ji as expected. 



ko/ki — s- limit - In the limit ki 3> kq 



become 



Atpo ~ ipo and Atpi ~ ipi, where e = kq/ki. Due 
to the strong screening in region 1, the potential under 
the surface wet by liquid 1 vanishes as ~ (Ti(eiKi)~^. 
The first iteration makes the approximation ip^^^ (r, 9) = 
/l^-* w 0. However, as a result of continuity of the poten- 
tial at 9 = a, and finite kq and co. 



shx( 



eoKQX cha;(7r — a)ch(7ra;/2) 



(12) 



Thus, cpoio^ > a) = -cr Vo"J^ 

in (p{P) arises from using ( JIS 
6* = a: 



0. The first nonzero term 
) in the jump condition at 



9./r(x, a) , = ^ (^^-^-^j^^ 

(13) 

Upon expanding i^^-c (^) about its dominant contribution 
at small ^ in the operator C^^ and performing the inte- 
gration over ^, 



2n 



60^0 tt> 



/■°° xT (n-Ux + x'))T(n+^{x- x')) e"''^' , 

Im / , , — r^— — — ^ : dx 

Jq ch(7rx72)chx'(7r - a)r(n + 1 - ix') 

_ (14) 
The contour integral ( [l4| ) can be performed exactly to 
find the two lowest order terms in e (for n = 0, fj,^^\e) = 
gir/(2ir-2a) ^(2) _ ^"j ^j-^g^^ must be Considered. 

Using (0), our final result is 



^(P;e/a^O) 
where 

To(a;£) = 



27r2 (o-o7o(a; e) + cr,„e) „ / ^ 2 



0(e^,e^), 
(15) 



^ _7r/(2ir-2Q) 



4(--«)^cos(^)r(||et) 



COS a 



(16) 



is plotted in Fig. || for various values of e. 

The kq/ki analysis fails for small a, when the fully 
screened approximation breaks down as the interface ap- 
proaches the cTi-charged solid substrate. The correction 
to the potential is independent of cri since it is screened 
out by a large ki. However, the jump condition across 
= a preserves an ei dependence. 

We have explicitly incorporated the geometric depen- 
dence of long-ranged electrostatic effects by deriving im- 
plicit equations for the contact angle valid for i?, Lh ~ 
{e/a)y/K-foi S> K^^- Our results show that for large 
mismatch in screening lengths, the potential at the con- 
tact point is proportional to the larger screening length 
but varies as an a-dependent power of the mis- 
match kq/ki. In the k » ki case, the zeroth order 
contribution to the potential can yield a nonnegligible 
effect: For k w 50nm, eq = 20, ei = 80, (Jq — cri w 0.1, 
(To, cTi ^ 0(1), the zeroth order correction to the surface 
energy varies from ~ 1.5 — OAksT/i^ as a varies from 




FIG. 3: The function Toia-e) (Eqn. (||)) that gives the ao- 
proportional term for the potential (Eqn. (|l5|)) in the £ = 
«;o/«:i limit. 



to TT. Since one ksT/i'^ ~ 8 dyne/cm, the implicit 
a dependence should be observable as long as the con- 
tact angle can be measured at a distance within 50nm 
of P. Such measurements are possible using AFM tech- 
niques pO[ |. Geometry-dependent electrostatics may also 
be an origin for the discrepancy between standard theory 
and measurements performed at low salt concentrations 
when screening length are long |8| . Even when the appar- 
ent contact angle cannot be measured within a screening 
length of P, our results may provide a basis by which 
to better understand dynamic phenomena. Since con- 
tact line pinning can occur at a length scale smaller than 
|P — P*|, the dynamics of the contact line may be better 
correlated with the true contact angle a rather than a* . 

Note that Si{a; A) and To(q;; e), when used in ( pi] ) al- 
low for the possibility of two solutions for a. The physical 
value for a will be determined by the minimum energy 
solution that needs to be determined by the full solution 
of the shape along P — P*. Therefore, in cases where two 
roots for a are possible, we expect the minimum energy 
solution to be the one closest to a* . Our results also im- 
ply that dramatic effects on wettability {a w 0, it) can be 
induced by small changes in the physical parameters. 

In our strictly 2D analyses, the only bounded solution 
when kq = ffQ ^ {e.g. air) is f{P; a) — 0. However, an 
asymptotic analysis for the large radius {R ^ Lh,K^^) 
limit is possible. The correction term to the Y-D equa- 
tion resulting from such an analysis is proportional to 
T (X (t^/(kJP) times an a-dependent factor, and defines 
the line tension 0, [2[ |l| arising from electrostatics. 
These electrostatic forces giving rise to line tension may 
be sufficiently long-ranged to be experimentally deter- 
mined through a measurement of a. Ionic strength may 
thus provide a controllable parameter with which to mea- 
sure the elusive, and controversial, line tension of wetting 
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